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We demonstrate a state reconstruction technique which provides either the Wigner function or 
the density matrix of a field mode and requires only avalanche photodetectors, without any phase 
or amplitude discrimination power. It represents an alternative to quantum homodyne tomography 
of simpler implementation. 
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I. INTRODUCTION 

The characterization of states and operations at the 
quantum level plays a leading role in the development of 
quantum technology. A state reconstruction technique 
is a method that provides the complete description of 
a physical system upon the measurements of an observ- 
able or a set of observables [l[ . An effective reconstruc- 
tion technique gives the maximum possible knowledge of 
the state, thus allowing one to make the best, at least 
the best probabilistic, predictions on the results of any 
measurement that may be performed on the system. At 
a first sight, there is an unavoidable tradeoff between 
the complexity of the detection scheme and the amount 
of extractable information, which can be used to recon- 
struct the quantum state [2[. Currently, the most ef- 
fective quantum state reconstruction technique for the 
radiation field is quantum homodyne tomography [H, 0| , 
which requires the measurement of a continuous set of 
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field quadrature and allows for the reliable reconstruc- 
tion of any Quantity expressible in terms of an expecta- 
tion value d,H, 0, B S G3 • A question arises on whether 
the tradeoff may be overcome by a suitable experimen- 
tal configuration or it corresponds to some fundamental 
limitations. Here we demonstrate that no specific dis- 
crimination power is required to the detector in either 
amplitude or phase, and that full state reconstruction is 
possible by a suitable processing of the data obtained 
with detectors revealing light in the simplest way, i.e. 
on/off detectors, such as single-photon avalanche photo- 
diodes. Of course, some form of phase and/or amplitude 
modulation is necessary, which, in our scheme, is im- 
posed to the field before the detection stage. In fact, 
our technique is built on the completeness of any set of 
displaced number states [H El El EI El El and the 
reliable maximum likelihood reconstruction of arbitrary 
photon- number distributions [l7| from on/off data. 

The paper is structured as follows. In Section |TT] we 
describe our reconstruction method, whereas in Section 
Mil the experimental setup used in the reconstruction is 
described in some details. Results are illustrated in Sec- 
tion llVl and the error analysis is reported in SectionfV] In 
Section [VI] we discuss few additional topics while Section 
IVIII closes the paper with some concluding remarks. 
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II. STATE RECONSTRUCTION BY ON/OFF 
MEASUREMENTS 



We start to describe our reconstruction technique by 
observing that the modulation of a given signal, de- 
scribed by the density matrix g, corresponds to the 
application of a coherent displacement (probe) g a — 
D(a)gD^ (a), a 6 C. In practice, it can be easily ob- 
tained by mixing the state under investigation with a 
known coherent reference jn a beam-splitter or a Mach- 
Zehnder interferometer [l8j |. Upon varying amplitude 
and phase of the coherent reference and/or the overall 
transmissivity of the interferometer, the modulation may 
be tuned in a relatively broad range of values. The main 
idea behind our method is simple: the photon distribu- 
tions of coherently modulated signals, i.e. the diagonal 
elements p n (ot) — {n\g a \n) of the density matrix g a , con- 
tain relevant information about the complete density ma- 
trix of the original signal g. Upon measuring or recon- 
structing the photon distribution p n (a) for different val- 
ues of the modulation one has enough information for full 
state reconstruction. By re-writing the above relation as 
Pn(a) = J2km D nk(a)gk m D mn {a), the off diagonal ma- 
trix elements may be recovered upon inversion by least 
square method, i.e. [l4[ 
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where TV^ is the number of modulating phases, n the 
truncation dimension of the Fock space, and F depends 
only on \a\ [14|. State reconstruction by the above for- 
mula requires, in principle, only phase modulation of the 
signal under investigation. Maximum likelihood methods 
and iterative procedures may be also used [19{ . On the 
other hand, the Wigner function may be reconstructed 
using its very definition in terms of displacement [20| 

W(a) = Tr[D{a) e rt(a) (-)"'»] = £(-)>„(a) . 



As a matter of fact, the measurement of the photon distri- 
bution is challenging as photo-detectors that can operate 
as photon counters are rather rare and affected either by 
a low quantum efficiency plj or require cryogenic con- 
ditions, thus impairing common use [22l. [23j. Therefore, 
a method with displacement but without photo-counting 
has been used so far only for states in the 0-1 subspace 
of the Fock space [l6j]. On the other hand, the exper- 
imental reconstructions of photon-number distributions 
for both continuous- wave and pulsed light beams is possi- 
ble using simple on/off single-photon avalanche photode- 
tectors. This requires the collection of the frequencies 
of the off events, P ,fc = J2™=o i 1 ~ Vk) n Pn at different 
quantum efficiencies of the detector, r\ k . The data are 
then used in a recursive maximum likelihood reconstruc- 
tion algorithm that yields the photon-number distribu- 



Pn +1 = Pn ^(A kn / A kj )(P , k /p , k [{pi}}), 
j 



k=l 



where A kn = (1 — T)k) and Po,fe is the probability of off 
events calculated from the reconstructed distribution at 
the zth iteration [24[. The effectiveness of the method 
has been demonstrated for single-mode [TtJ and multi- 
mode fields [25l ] . and also applied to improve quantum 
key distribution [26| . 

Since the implementation of the modulation is rel- 
atively easy, we have thus a reconstruction technique 
which provides the quantum state of radiation modes and 
requires only avalanche detectors, without any phase or 
amplitude discrimination power. Here, we develop the 
idea into a proper reconstruction technique and demon- 
strate the reconstruction of the Wigner function [27[ and 
the density matrix for different states of the optical field. 



III. EXPERIMENTAL SETUPS 

We have performed two experiments for the recon- 
struction of the Wigner function and the density matrix 
respectively. In Fig. [T] we sketch the corresponding ex- 
perimental setups: the upper panel for the measurement 
of the Wigner function and the lower panel for the den- 
sity matrix (lower panel) . The first set of measurements 
was performed on ps-pulsed light fields at 523 nm wave- 
length. The light source was the second-harmonic output 
of a Nd:YLF mode-locked laser amplified at 5 kHz (High 
Q Laser Production) delivering pulses of ~ 5.4 ps dura- 
tion. 



Spatial filter 




Spatial filter 



FIG. 1: (Color online) Schematic diagram of the two exper- 
imental setups. Upper panel: Nd:YLF, pulsed laser source; 
P, polarizer; HPD, hybrid photodetector; SGI, synchronous 
gated integrator; ADC, analog-to-digital converter. Lower 
panel: He-Ne, cw laser source; IF, interference filter; APD, 
single-photon avalanche photodiode. F, neutral-density fil- 
ter; BS, beam splitter; Pz, piezoelectric movement; R, ro- 
tating ground glass plate. Components in dotted boxes are 
inserted/activated when necessary. 

The field, spatially filtered by means of a confocal tele- 
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scope, was split into two parts to give both signal and 
probe. The photon-number distribution of the probe was 
kept Poissonian, while the coherent photon-number dis- 
tribution of the signal field was modified in order to get 
suitable states of light. In particular, we have generated 
two phase-insensitive classical states, namely a phase- 
averaged coherent state and a single-mode thermal state. 
The first one was obtained by changing the relative phase 
between signal and probe fields at a frequency of ~ 300 
Hz by means of a piezoelectric movement (Pz in the up- 
per panel of Fig. [I}, covering 1.28 /zm span. On the other 
hand, to get the single-mode thermal state we inserted 
an Arecchi's rotating ground glass disk on the pathway 
of the signal field followed by a pin-hole that selected a 
single coherence area. 

Signal and probe fields were mixed in an unpolariz- 
ing cube beam-splitter (BS) and a portion of the exiting 
field was sent, through a multimode optical fiber (600 
/im core-diameter), to a hybrid photodetector module 
(HPD, mod. H8236-40, Hamamatsu, maximum quantum 
efficiency ?7hpd=0.4 at 550 nm). Although the detector 
is endowed with partial photon resolving capability, we 
used it as a simple on/off detector by setting a threshold 
at the value corresponding to zero detected photons. Its 
output current pulses were suitably gate-integrated by 
a SR250 module (Stanford Research Systems, CA) and 
sampled to produce a voltage, which was digitized and 
recorded at each shot. In order to modulate the ampli- 
tude of the probe field, a variable neutral density filter 
was placed on its pathway. The maximum overall detec- 
tion efficiency, calculated by including the losses of the 
collection optics, was r] max = 0.29. We used a polarizer 
put in front of the fiber to vary the quantum efficiency 
of the detection chain from rj max to 0. 



IV. RESULTS 

Here we illustrate the reconstruction obtained for 
the Wigner function and the density matrix of phase- 
averaged coherent states and thermal states, as obtained 
by our method after recording the on/off statistics of 
amplitude- and/or phase-modulated signals. 

In Fig. [5] we report the reconstructed Wigner functions 
for a phase-averaged coherent state with real amplitude 
z ~ 2.1 and a thermal state with average number of pho- 
tons nth — 2.4. The Wigner function of the vacuum state 
is also reported for comparison. As it is apparent from 
the plots all the relevant features of the Wigner functions 
are well recovered, including oscillations and the broad- 
ening due to thermal noise. In this case raw data are 
frequencies of the off event (collected over 30000 laser 
shots) as a function of detector efficiency, taken at differ- 
ent amplitudes of the modulating field, whereas the in- 
termediate step corresponds to the reconstruction of the 
p n {a). The insets of Fig. [2] display p n {a) for the phase- 
averaged coherent state at two values of the modulation. 
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FIG. 2: (Color online) State reconstruction by amplitude- 
modulation and on/off measurements. In the main plot we 
report the off frequencies as a function of the quantum effi- 
ciency as obtained when the signal under investigation is a 
phase-averaged coherent state and for different values of the 
probe intensity |a| 2 . The two insets show the reconstructed 
photon distributions for \a\ 2 — 5.02 and \a\ 2 = 10.69, cor- 
responding to the off distributions indicated by the arrows. 
The vertical black bars denote the mean value of the photon 
number for the two distributions ((o T a) = 11.3, upper, and 
(aJa) = 18.0, lower). In the lower left plot we report the 
corresponding reconstructed Wigner function. In the lower 
right plot we report the Wigner functions for signals in (blue) 
thermal state and (yellow, with sharper peak) vacuum . 



The second set of measurements has been performed 
to achieve state reconstruction with phase modulation. 
Here the light source was a He-Ne laser beam attenu- 
ated to single-photon regime by neutral density filters. 
The spatial profile of the beam was purified from non- 
Gaussian components by a spatial filter. Beyond a beam- 
splitter, part of the beam was addressed to a control 
detector in order to monitor the laser amplitude fluc- 
tuations, while the remaining part was sent to a Mach- 
Zehnder interferometer. A piezo-movement system al- 
lowed changing the phase between the "short" and "long" 
paths by driving the position of the reflecting prism with 
nanometric resolution and high stability. The beam on 
the short arm represented the probe, while the beam 
on the long arm was the state to be reconstructed. In 
the first acquisition the signal was the coherent state 
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FIG. 3: (Color online) State reconstruction by phase- 
modulation and on/off measurements. In the upper plot we 
report the off frequencies as a function of the quantum ef- 
ficiency as obtained when the signal under investigation is 
a coherent state and for different phase-shifts. The two in- 
sets show the reconstructed photon distributions for the two 
phase-modulated versions of the signal corresponding to max- 
imum and minimum visibility at the output of the Mach- 
Zehnder interferometer. The vertical black bars denote the 
mean value of the photon number for the two distributions, 
(a^a) = 3.5 and (a^a) — 2.9. In the lower left plot (left) we 
report the corresponding reconstructed density matrix in the 
Fock representation (diagonal and subdiagonal elements). In 
the lower right plot we report the density matrix for the signal 
excited in a thermal state. 



off event as a function of the detector efficiency, taken at 
different phase modulations, (f>, whereas the intermediate 
step corresponds to the reconstruction of the photon dis- 
tribution for the phase-modulated signals. In our exper- 



iments we used Na, — 12 and 



= 0.01 for the coherent 



state and 



= 1.77 for the thermal state. The use of a 



larger N$ would allow the reliable reconstruction of far 
off-diagonal elements, which is not possible in the present 
configuration. In the insets of Fig. [3] we report the recon- 
structed distributions at the minimum and maximum of 
the interference fringes. 



ERROR ANALYSIS 



The evaluation of uncertainties on the reconstructed 
states involves the contributions of experimental fluctu- 
ations of on/off frequencies as well as the statistical fluc- 
tuations connected with photon-number reconstruction. 
It has been argued [H, [29| that fluctuations involved in 
the reconstruction of the photon distribution may gen- 
erally result in substantial limitations in the information 
obtainable on the quantum state, e.g. in the case of mul- 
tipeaked distributions [3(| ■ For our purposes this implies 
that neither large displacement amplitudes may be em- 
ployed, nor states with large field and/or energy may be 
reliably reconstructed, although the mean values of the 
fields measured here are definitely non-negligible. On the 
other hand, for the relevant regime of weak field or low 
energy, observables characterizing the quantum state can 
be safely evaluated from experimental data, including e.g. 
the parity operator used to reconstruct the Wigner func- 
tion in the phase-space. In our experiments, the errors 
on the reconstructed Wigner function are of the order of 
the size of the symbols in Fig. [3] whereas the absolute 







exp 



Hi 



on the reconstruction of the 



density matrix in the Fock basis are reported in Fig. 0J 



itself while in the second acquisition it was a pseudo- 
thermal state generated as described above. The detec- 
tor, a Perkin-Elmer Single Photon Avalanche Photodi- 
ode (SPCM-AQR) with quantum efficiency r\ max = 0.67, 
was gated by a 20 ns time window at a repetition rate 
of 200 kHz. To obtain a reasonable statistics, a single 
run consisted of 5 repetitions of acquisitions lasting 4 s 
each. In the bottom part of Fig. [3] we report the recon- 
structed density matrix in the Fock representation (diag- 
onal and subdiagonal) for a coherent state with real am- 
plitude z ~ 1.8 and a thermal state with average number 
of photons equal to n t h — 1.4. 

As it is apparent from the plots the off-diagonal ele- 
ments are correctly reproduced in both cases despite the 
limited visibility. Here the raw data are frequencies of the 




FIG. 4: (Color online) Absolute difference A nm = |£Vm — f?rim| 
between reconstructed and theoretical values of the density 
matrix elements for the coherent (left) and thermal (right) 
states used in our experiments. 
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VI. DISCUSSIONS 

We have so far reconstructed the Wigner function and 
the density matrix for coherent and thermal states. The 
extension to highly nonclassical states does not require 
qualitative changes in the setups. The only difference 
stays in the displacement, which should be obtained with 
high transmissivity beam splitter in order to avoid mixing 
of the signal [l8| . 

As our method involves a beam splitter where the sig- 
nal interferes with a reference state in order to obtain 
the displacement, we have optimized the effectiveness of 
mode matching and of the overall scheme by standard 
visibility test. We note that in this point our technique 
is similar to quantum homodyne tomography (QHT), 
where the signal is amplified by the mixing at a beam 
splitter with a strong local oscillator. The main differ- 
ence with standard QHT, however, is the spectral do- 
main of the measurement, which for QHT is confined to 
the sideband chosen for the measurement, while it is not 
in our case. The use of pulsed temporal homodyning 
[HI H, [H, HI SH would remove this limitation of QHT. 
However, this technique is still challenging from the ex- 
perimental point of view and thus of limited use. The 
effect of photodetectors efficiency should be also taken 
into account. This is a crucial issue for QHT, which may 
even prevent effective reconstruction (3(| . For the present 
on/off reconstrucion, it does not dramatically affect the 
accuracy 24J. Notice also that any uncertainty in the 
nominal efficiency rjmax of the involved photodetectors 



does not substantially affect the reconstruction [2J|. 

We stress that our method is especially suited for low 
excited states, as it does not involve intense fields and 
delicate balancing to reveal the relevant quantum fluctu- 
ations. 



VII. CONCLUSIONS 

In conclusion, we have demonstrated a state recon- 
struction technique providing Wigner function and den- 
sity matrix of a field mode starting from on/off photode- 
tection of amplitude- and/or phase-modulated versions 
of the signal under investigation. Our scheme is lit- 
tle demanding as to the detectors, with the amplitude 
and phase control required for full state characterization 
transferred to the optical setup, and appears to be re- 
liable and simple especially for states with small num- 
ber of photons. We foresee a possible widespread use in 
emerging quantum technologies like quantum informa- 
tion, metrology and lithography. 
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